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Abstract 

Fractional action-like variational problems have recently gained 
importance in studying dynamics of nonconservative systems. In this 
note we address multi-dimensional fractional action-like problems of 
the calculus of variations. 
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1 Introduction 



Fractional calculus (FC) represents a powerful tool in applied mathematics to 
study a myriad of problems from different fields of science and engineering, 
with many break-through results found in mathematical physics, finance, 
hydrology, biophysics, thermodynamics, control theory, statistical mechanics, 
astrophysics, cosmology and bioengineering [2], [20l [211 1211 [251 [SHI [271 [301 
[3T1 l32j . Although various fields of application of fractional derivatives and 
integrals are already well established, some others have just started. An 
example of the later is the study of fractional problems of the calculus of 
variations and respective Euler-Lagrange type equations [H [H [SI [HI [II]- 

During last years, fractional variational principles (FVP) were proposed 
and several applications given. Different methods were used to obtain the 
fractional Euler-Lagrange equations and the corresponding Hamiltonian canon- 
ical equations [Il[3l[2[5l[6l[Tl[22l[23l[28l[29l[33]. The existence of several 
different FVP and the need for still more elaboration on the subject, both 
for classical and quantized systems, is partially explained by the nonlocal 
nature of the fractional differential operators used to describe the dynamics, 
and the form of the corresponding adjoint operators. Another reason is the 
existence of many different fractional integral operators, including the ones of 
Griinwald-Letnikov, Caputo, Riesz, and Riemann-Liouville. The Riemann- 
Liouville operator is one of the most frequently used when fractional integra- 
tion is performed |2T1 [30] . 

Recently, the first author has proposed a new one-dimensional (ID) ap- 
proach, known as the fractional action- like variational approach (FALVA), 
in order to model nonconservative dynamical systems [3 [S]. In FALVA the 
fractional time integral introduces (only) one real parameter a, and the de- 
rived Euler-Lagrange equations are simpler and similar to the standard ones. 
The novelty in the Euler-Lagrange equations is the presence of a fractional 
generalized external force acting on the system [71 [H [9]. In particular, no 
fractional derivatives appear in the derived equations. The conjugate mo- 
mentum, the Hamiltonian and Hamilton's equations are shown to depend 
on the fractional order of integration a. Constants of motion for fractional 
action-like variational problems were discussed in [151 [16] (see also [T71 [18]): 
FALVA problems with higher-order derivatives are studied in [T6l [19] : and 
some encouraging results obtained and discussed in [TOl [Til [I2l [111 E] . Here 
we are interested to generalize FALVA for multiple fractional action-like in- 
tegrals of the calculus of variations (multi-dimensional fractional action-like 
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variational problems). 

The text is organized as follows: in Sect. [2] we review the basic concepts 
of ID-FALVA; the new extensions are found in Sect. [31 Since various applica- 
tions of FC are based on replacing the time derivative in dynamical equations 
with a fractional derivative, in §3.11 we introduce the double-weighted EVP 
with the recent fractional derivatives of J. Cresson [B]. In §3.21 the FVP 
is given in general form, for the arbitrary A^-dimensional case. Sect. S] is 
dedicated to conclusions and future perspectives. 

A Note on the notation used 

In order to be clear when f{t) stands for a function or the value of the 
function at a point t, we denote the function hy t —>■ f{t) and the value of 
the function at t by f{t). In the notation t f{t), t is a dummy variable. 
For instance, in the beginning of §2.21 we write {q,q,T) —>■ L{q,q,T), which 
represents a function of three variables. Exactly the same function can be 
written, for example, as (a, b, c) L{a, b, c) {q, q, r or a, b, c are here dummy 
variables). However, the dummy variables we choose in the paper are used 
consistently to fix the notation for the partial derivatives. For example, if 
we define the function as (a, c) L(a, 6, c), then we write ^ to denote 
the partial derivative of function L with respect to the first argument; if we 
define it as (g, g, r) L{q, q, r), then the partial derivative of L with respect 
to the first argument is denoted by 

2 Brief Overview of ID-FALVA 

In this section we summarize the one-dimensional FALVA. The review con- 
sists of two cases: absence of fractional derivatives ( §2.ip : and presence of 
Riemann-Liouville fractional derivatives in the sense of Cresson [6] ( §2.2p . 
The reader is referred respectively to [71, fS] and [H] for more details. 

2.1 Absence of fractional derivatives 

Consider a smooth n-dimensional manifold M (configuration space) and de- 
note by L : TM x M x M — M the smooth Lagrangian function. For any 
smooth path q : [a,b] ^ M satisfying fixed boundary conditions q{a) = qa 
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and q{b) = Qf,, we define the fractional action integral by 



t 



S''[q]{t) = Y^J mr),q{r),r){t-rr-'dr, (1) 

a 

where T is the Euler gamma function, g = < a < 1, r G (a, t) is 
the intrinsic time and t G [a, b] is the observer time. For a discussion of 
the importance to consider a multi-time formalism, we refer the reader to 
[311 [35]. The Lagrangian L{q, q, r) is weighted by (t — r)°~Yr(a), and one 



can write the smooth action integral 5*" [q\ (t) as / L{q{r), qi^r), T)dgt{T), with 

a 

an appropriate time smeared measure dgt{r) on the time interval [a,t], or 
as the Riemann-Liouville operator applied to the Lagrangian L{q' [t) , q{t) , t) 
[3, [H]. Here we just note that when a ^ 1 the functional S°' [q] (6) tends to 
the classical action integral of the calculus of variations: 





]im S^[q]{b)= [ L{q{r),q{r),T)dT. 



Functionals of type ([Tj) appear naturally in mathematical economy, where 
they are used for describing discounting economical dynamics; and in the 
theory of dynamical systems, for describing nonlinear dissipative structures. 

The Euler-Lagrange equations associated with the fractional action inte- 
gral ([T|) take the following form: for all t G [a, b], 

dL_d_ /dL\ _ 1-a dL . _ ^ 
% dr \dqij t-r dqt' 

where the partial derivatives of L are evaluated at (g(r), g(r), r), r G (a,t). 
The fractional term ^^f^ has a meaning in classical, quantum or relativistic 
settings, being responsible for adding a time-dependent damping coefficient 
into the dynamical equations. If we denote by i? = {1 — a)L/{t — t) the frac- 
tional Rayleigh dissipation function and by E the Euler-Lagrange operator, 
i.e. E = — then equations ([2|) take the form E{L) = Extremals 
are defined as solutions of the fractional Euler-Lagrange equations (FELE) 
EiL) = f . 
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2.2 Presence of fractional derivatives 



For any smooth path g : [a, 6] — > M satisfying fixed boundary conditions 
g(a) = Qa and q{b) = qt,, let the fractional functional associated to (g, q, r) — *• 
L(g, g, r) be now defined by 



L(D:;'^g(r),g(r),r)(t-r)"-Mr, 



(3) 



where [a, t] C [a, 6] C M, and D"''^ denotes the Riemann-Liouville derivative 
of order (a,/9) as defined by J. Cresson in [6j: for all a, t G M, a < t, the 
fractional derivative operator of order (a, /9), < a, < 1, is given by 



«7 

y 



(4) 



with 7 G C, i = V— 1, D"^ and Df_ the usual left and right Riemann- 
Liouville fractional derivatives of order < a,/? < 1 [30] , 



d 



rri-a) de 



f{r) {e-ry^dr, 



1 



d 



T{i-p) \ de 



f(r){r-9y^dT\ , 



9 G [a, t]. For more details we refer the reader to [6]. Here we just note that 
for 7 = i one has D"'^ = — -Df_; for 7 = — i, Z?"'^ = i^^i ; for a — 1 and 



/3 — i> 1 one obtains 1^°'^ = ^ and ([3]) is reduced to the classical functional of 
the calculus of variations. In 
to ([3]) is given by 

dL 



it is proved that the FELE correspondent 



dq 



-D 



I3,a 



fdL 



- 



1 — a dL 
t — T dq 



(5) 



V r G (a, t) , where the partial derivatives of L are evaluated at {D"'^q{r) ,q{r), r) 

In the next section we generalize ([5]) for FALVA problems with multiple 
integrals. 
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3 Main Results 



More generally, we discuss here the multi-dimensional fractional variational 
problem or A^D-FALVA (A^-dimensional FALVA). We first start by double 
weighted fractional integrals, i.e. by the 2D-FALVA. 



3.1 Double- Weighted Fractional Variational Principles 
(2D-FALVA) 

Similar to Sect. [2], we denote by M a smooth n-dimensional manifold; the ad- 
missible paths are smooth functions q : Q cM."^ ^ M satisfying fixed Dirich- 
let boundary conditions on dfl; the Lagrangian function {qj.,qy,q,x,y) — >■ 
L{qx,qy,q,x,y) is supposed to be sufficiently smooth with respect to all its 
arguments; a and p are two real parameters taking values on the interval 
(0,1). 

Definition 1. The 2D-FALVA action integral is defined by 

Y^J^^fT^^ 11^^^ L{D^iQ, D^^l^Q, <1, ^, y) (e - xT"' (A - yf-' dxdy , (6) 

where and A are the observer times, A) G ^2; x and y are the intrinsic 
times, {x,y) E ^1{^,X) C f2, x 7^ ^ and y ^ X; q = q{x,y)] D^i" and 
D^:^ are the fractional derivative operators (jl]), respectively of order {S,a) 
with respect to x and of order (x, (3) with respect to y. We denote ([6]) by 
S^'P'^'^ [q] (e,A). 

Remark 1. The classical multi- variable variational calculus has some limita- 
tions which a multi-time variational calculus, like the one we are proposing 
here, can successfully overcome. For a discussion of the importance to con- 
sider a multi-time formalism we refer the reader to the works of C. Udriste 
and his collaborators [3ll [35] . 



The primary objective is to find functions q = q{x, y) that make the 
fractional action S""'^'^'-^ [q] (^, A) stationary for every (^, A) G Q. 

Theorem 1. Given a smooth Lagrangian {qx,Qy,<l,x,y) —>■ L{qx,qy,q,x,y), 
'i'f Q = q{x^ y) makes the fractional action S""'^'^'^ [q] (^, A) defined by ([6]) 
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stationary for every A) e ft, then the following double-weighted Euler- 
Lagrange equation holds for every {x, y) e A); 

^5,a ( dL{D^iq,D^^:^q,q,x,y) \ . ( dLjP'^^q, D^^:^q,q,x,y ) 



( dL{D-f^q, D^^i^q, q,x,y)\ 1-^ f dL{D'^f^q, D^^i^q, q, x, y) 



i-x \ dq^ ) X-y \ dqy 

.y) 

0. (7) 



dL{D^f^q,D^:;^q,q,x,y) 



dq 

Proof. Let g be a stationary solution, /?. ^ 1 a small real parameter, and 
w{x,y) an admissible variation, i.e. an arbitrary smooth function satisfying 
w{x,y) — for all {x,y) e dO, so that q + hw satisfies the given Dirichlet 
boundary conditions for all h. The fractional action S^'^'^'^ [q + hw] can be 
written as 

— / / L{D'^!iq + hD^!iw, D^^i^q + q + hw, x, y) 

X (e - xy-' (A - yf-' dxdy , 
and the stationary condition -^iS^'^'^'^ [q + hw]\^_^ = gives 

1 ff f dL dL dL 

w— + D":iw— + D^:^w— 



r(a)r(/?)yya(j,,) V dq ^-'-^ dq, ^'^ dqy, 

xi^-xr-'{X-yf-'dxdy = 0, (8) 

where the partial derivatives of function {qx,qy,q,x,y) — > L{qx,qy,q,x,y) 
are evaluated at {D^:^q{x,y),D^:^q{x,y),q{x,y),x,y). Using integration by 
parts and Green's theorem, we know that 



f2(e,A) 



fdP dP \ - - 



for any smooth functions Gi and G2, where 

r(l + «)A = A^-(A-yf 



We conclude from that 



r(a)r(/3)7^(^,,) 



dqx) 



) (e - xT-' (A - yf-' + (1 - /?) ( ^ ) - x)"-^ (A - y) 



\/3-l 



\/3-2 



(ixc/?/ = 



and then, because of the arbitrariness of w{x, y) inside A), it follows ([7]), 
which is the 2D-FELE. □ 

We expect that fractional variational problems involving multiple inte- 
grals may have important consequences in mechanical problems involving 
dissipative systems with infinitely many degrees of freedom. 



3.2 N- Weighted Fractional Variational Principles (ND- 
FALVA) 

All the arguments of §3.11 can be repeated, mutatis mutandis, to the A^- 
dimensional situation when the admissible paths are smooth functions q : 
VldM.^ ^ M satisfying given Dirichlet boundary conditions on dVt. 

Definition 2. Consider a smooth manifold M and let {qx^ , • • • , qx^ A^^i^ ■ ■ ■ ^ ^n) 
L{qx^, . . . , qxj^, q,xi, . . . ,xn) be a sufficiently smooth Lagrangian function. 
The AD-FALVA functional is defined by 

N 

Sf[q]iO = ir^ / ■■■ / L{V^/qix),q{x),x)ll{^.-x.r-'dx, 

n r (a,) ^=1 

i=l 

where x = {xi, . . . ,xn) is the intrinsic time vector, ^ = ■ ■ ■ ,^n) £ ^ 
the observer time vector, x G C fl with Xi ^ {i = 1, . . . ,N), dx = 
dxi ■ ■ ■ dx^, a = (ai, . . . , Oat), S = {6i, . . . , 6^), < ctj < 1 (i = 1,...A^), 
and = iD-}/^\...,D^!^,'^^). 
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Theorem 2. The N -dimensional Euler- Lagrange equation associated to the 
fractional functional S"'^ [q\ {^), ^ e is given by 



N 

E 

i=l 



D 



dL 
dqxi 



1 — a,- / dL 



dL 

dq 



= 0, 



where all partial derivatives of L are evaluated at (V^'^q{x),q{x), x) , x e 



4 Conclusions and Further Work 

In this work wc introduce the muhi-dimensional FALVA setting and derive 
the corresponding muhi-dimensional fractional Euler-Lagrange equations. 
Obtained Euler-Lagrange equations are enough complicated, and one expects 
solutions to be found using numerical techniques. A work in this direction is 
in progress. We expect that fractional variational problems involving multi- 
ple integrals will have important consequences in mechanical problems and 
optimal control theory involving dissipative systems with infinitely many de- 
grees of freedom. A geometric formulation of the field equations for the 
fractional action-like variational formalism, in terms of multi-vector fields on 
tangent bundles, is under investigation. 

In our paper, as well as in all previous works on fractional Euler-Lagrange 
equations we are aware of, it is assumed that at least one stationary point 
for the fractional functional exist. Euler-Lagrange equations are valid under 
this assumption. The question of obtaining conditions on the Lagrangian L 
assuring the existence of stationary trajectories is, to the best of our knowl- 
edge, a completely open question in the fractional setting. This is a pertinent 
question because even very simple Lagrangians, e.g. L = D^'^q, fail to satisfy 
the hypotheses under which the Euler-Lagrange equations are valid. 
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